Nanomechanical read-out of a single spin 
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c/2 Abstract 

a 

The spin of a single electron in a suspended carbon nanotube can be read out by using its 

coupling to the nano-mechanical motion of the nanotube. To show this, we consider a single 

electron confined within a quantum dot formed by the suspended carbon nanotube. The spin- 

orbit interaction induces a coupling between the spin and one of the bending modes of the 

suspended part of the nanotube. We calculate the response of the system to pulsed external 

J> driving of the mechanical motion using a Jaynes-Cummings model. To account for resonator 
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damping, we solve a quantum master equation, with parameters comparable to those used in 
recent experiments, and show how information of the spin state of the system can be acquired 
by measuring its mechanical motion. The latter can be detected by observing the current 
through a nearby charge detector. 

In the past few years, advanced manufacturing techniques have spawned new interest in nano- 
mechanical devices. These are of interest for fundamental research as well for possible applica- 
tions. In various experimental setups the cooling of a nanomechanical resonator to its ground state 
has been achieved.^' It is now possible to study quantum features such as the zero-point motion of 
objects much larger than the atomic scale. 4 As mechanical motion can be coupled via a wide range 
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of forces, nanomechanical systems have been proposed for various applications such as mass,^ 
force, 7 and motion^ sensing. This versatility also makes hybrid systems of mechanical devices 
coupled to other systems possible candidates for various applications in quantum information and 
communications. For example, nanomechanical resonators have been proposed as qubits, 9 optical 
delay lines,^ quantum data buses, 11 and quantum routers 12 among others. 

One of the challenges in building sensitive nano-electromechanical systems is to control ther- 
mal fluctuations. Only if the device operates at energies significantly higher than the thermal 
energy, one can expect to observe quantum mechanical behavior in equilibrium. In this context, 
carbon is a very promising building material because the resulting structures are very light and 
stiff which leads to high resonance frequencies. Resonators made of suspended carbon nanotubes 
(CNTs) are of special interest. They can be produced almost free of defects and with radii on 
the scale of ~ 1 nm and lengths up to 1 jUm and in addition their bending mode is easily excited. 
Experimentally it has been demonstrated that frequencies of more than 4 GHz and quality factors 
of more than 100,000 are achievable J 13 l 14 l 

Coupling a localized electron spin to a nanomechanical resonator is of particular interest be- 
cause such a system may be used as a quantum memory due to the relatively long spin lifetimes.^ 
Carbon based materials with their weak spin-orbit coupling due to the relatively low atomic mass 
of carbon and the presence of only few nuclear spins is a promising material.^ Recently the read- 
out of a single spin of the NV center in diamond has been used to read out the oscillatory motion 
of the magnetized tip of an atomic force microscope.^^ However, the well-defined local field 
gradient required for this setup is quite challenging. CNTs provide an entirely different arena for 
the coupling of single spins to nanomechanical resonators and they can potentially be integrated 
into scalable (2D) semiconductor structures. 

Here, we propose a method for an all-electrical readout of the state of a single electronic spin. 
We assume that the spin is confined to an oscillator consisting of a doubly clamped, suspended 
CNT While previous experiments^ 13 ! 14 ! make use of a charge transport measurement throught the 
CNT, here we focus on the measurement of the amplitude of the vibrating CNT via a close-by 
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Figure 1 : Schematic of the proposed setup of the carbon nanotube (CNT) resonator with a quantum 
dot containing a single electron. The gates on both sides of the CNT can be electrostatically 
adjusted to control the number of confined electrons. The backgate is used to tune the resonance 
frequency of the CNT into resonance with the splitting of the Zeeman sublevels of the electron 
which make up the qubit. The charge sensing device in the vicinity to the vibrating nanotube 
serves as a detection device for changes in the amplitude due to the coupling of the spin to the 
vibrational motion. It can be a quantum point contact (QPC) or a quantum dot (QD) which needs 
to be operated in a regime where the current depends non-linearly on the gate voltage. A current 
meter is used to measure a nonlinear current / which flows through the device. 

charge sensor, which could be a quantum point contact (QPC) or a quantum dot (QD). Those 
devices have been shown to be sensitive to charge fluctuations of far less than one elementary 
charge.^^The coupling of the spin to the vibrational motion allows the read-out of the spin state 
via measurements of the amplitude of the mechanical vibration. 23 This scheme does not require 
optical access to the probe or time-dependent magnetic fields which makes it potentially scalable. 

The setup under study is shown in Figure ([I]). It consists of a CNT suspended over a trench 
with contacts on both sides and electrostatic gates responsible for the electron confinement. Since 
the nanotube carries a net charge, a constant (dc) voltage applied to a back gate can be used to tune 
the mechanical resonance frequency. Furthermore, an ac voltage applied to an external antenna or 
the back gate can be used to drive the motion of the charged resonator J 14 * 23 l 

In order to obtain a well defined qubit, the fourfold degeneracy due to the valley and spin 
degree of freedom of the CNT has to be lifted (note that for the localized electron, the sublattice 
degree of freedom is frozen out). A magnetic field B parallel to the CNT axis in combination with 
the spin-orbit coupling serves this purpose. Around a field of B* « Aso/^b two orthogonal spin 
states within the same valley are split by the Zeeman energy hco q = /is(5 — B*). Here Aso refers 
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to the intrinsic spin-orbit coupling strength; the intervalley coupling A KK / is assumed to be much 
smaller. In the following we assume A$o = 370 /leV and A KK / = 65 /leV as realistic parameters.^ 
The corresponding states in the opposite valleys are energetically well seperated. Thus we can 
treat the system as an effective two level system. 

For the detection of the nano-mechanical oscillation we propose to use that the position x(t) 
of the oscillating charged CNT modulates the current through the charge sensor in a (generally) 
nonlinear way, i.e. I(t) = Iq + I\x(t) + I2X (t) + 0(x 3 ). This is the case, for example, when a 
QD is operated at the maximum of a Coulomb blockade peak. This situation will be assumed for 
following discussion. Furthermore, we assume that the frequency co p of the oscillator is much 
larger than the tunneling rate through the charge sensing device. In this case, in lowest order 
the latter only probes a time averaged squared displacement X 2 = }- J Q Z dtX 2 (t) where X 2 (t) = 
(x 2 (t)) = Tr(x 2 p(t)) with the oscillator density matrix p. In the expansion of I(t), only even-order 
terms in x appear because x n = for n odd. Higher order terms n > 2 are neglected because they 
contribute only weakly to the current for small displacements x(t) . Our goal in this paper is to 
calculate the time-averaged current 7 through the charge sensing device as a function of the spin 
state of the electron. 

In the following, we restrict our considerations to one polarization of the bending mode of 
the CNT. The generalization to two modes is straightforward. When quantized, the resonator dis- 
placement x can be written as x = {a + a*) , where a and are phonon creation and annihilation 
operators and Iq is the zero-point motion amplitude of oscillator. In principle there are two ways in 
which the flexural phonons can couple to the spin. At large phonon energies the usual deformation 
potential is dominant as it is proportional to q 2 where q is the phonon wave number. For lower 
energies however the spin-orbit mediated deflection coupling q to flexural phonons dominates.^ 
The resulting coupling strength g will be proportional to both the spin-orbit coupling A50 and the 
zero-point amplitude IqJ^ 
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The system is described by the Jaynes-Cummings Hamiltonian 

H = H + H d (1) 

with 

Ho = —^-o z + hg(a(7 + +a^(7-) +ha) p a f a, Hj = hX(a + a'), (2) 

where o z is a Pauli matrix acting on the spin qubit while a and a* are the creation and annihilation 
operators for the phonons. The qubit and oscillator frequencies in the rotating frame, (b q = (O q — (O 
and cbp — COp — ft), are given as detunings from the driving frequency ft). The driving strength A is 
assumed to be weak, i.e., X <^ C0 p . In the Jaynes-Cummings model a rotating-wave approximation 
is incorporated which is valid as long as the driving frequency ft) is comparable to C0 q and Q) p . 

In addition to the unitary evolution we include the damping of the CNT with a rate T which can 
occur on the same timescale as the read-out via the charge detection device. The spontaneous qubit 
relaxation y=\/T\, where T\ denotes the spin relaxation time, is neglected because of the very low 
density of other phonon modes in the vicinity of the bending mode at frequency (O p . Previously^ 
we have shown that with state-of-the-art experimenental techniques the strong-coupling regime, 
i.e. g 3> r, 7, is within reach. The non-unitary dynamics are described by the quantum master 
equation for the time evolution of the density matrix p 

+ (n B + l)r(apa t -^{a 1 a,p}) (3) 
+ ngr(a^pa — ~{aa\p}). 

Here ng = \j[e h<S) pl kBT — 1) refers to the Bose-Einstein occupation factor of the phonon bath at 
temperature T. The Lindblad terms T correspond to emission (absorption) of a phonon to (from) 
the phonon bath. 

While solutions of the quantum master equation (??) cannot be given in closed form, it is 
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worthwhile studying the eigenstates of Hq first, 



|V„,+> =cos-| t«-l) + sin-||n), (4) 
|W) = - sm ylt 1/ + COS — | 

for n > l and the special case of the ground state n = 0: | yb) = |0 4)- Here, we use the notation 
\no) for eigenstates of Hq withg = 0, i.e., a^a\no) =n\na) and o z \no) = o\no) and o =t,|= ±1. 
The mixing angle a is defined by tana = ? 8v ^ ■ In the resonant case a — > 7t/2. The eigenenergies 
are 

E n ,± = hco p (n - ^)±^(^y^) 2 +^ (5) 

for n > 1 and Eq = —h^. For n = 0, we have | \. 0} as the ground state. The energy splitting 
between adjacent eigenstates (see Figure (|2^)) in the resonant case cb p = (D q is E n +i_± —E n> ± = 
± (h(b p + hg(y/n+ 1 — a/h)) and E n>+ — E n _ = 2hgy/n. 

Since we are interested in the read-out of the electron spin, we assume an initial state in which 
the latter is in one of the two o z eigenstates. In the case of T = and an empty QD the osciallator 
is in its ground state |0) with the zero-point amplitude Iq. Directly after loading with an elec- 
tron the state of the system is |0 t) or |0 \). At temperatures T > 0, the distribution of phonons 
obeys Bose-Einstein statistics and we obtain y¥(a))r = |c) <8> \Y^=Q e ~ nhCOp bT \ n ) as a initial 
state, where o =t,4 and Z = YZ=Q e ' lh(i> ^ kBT l% me partition function. Experimentally, the state 
preparation could be performed using techniques that have been established in conventional semi- 
conductor QDs.^ Raising one of the gates well above the two Zeeman-split qubit levels allows 
an electron to hop onto the CNT Provided geiMs^ ^> kgT.g, where g e i is the electron g-factor, /ig 
the Bohr magneton, and B the magnetic field, the probability for finding either Zeeman sublevel 
to be populated is 1/2. A subsequent measurement will then reveal which state was prepared. If 
the measurement is delayed by a time T, then the population of the higher-energy spin state will 
decay as ~ exp(— t/7i)J^To demonstrate the coupling between spin and phonon and to study the 
read-out this random filling method is sufficient. For possible application in quantum computing, 
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Figure 2: (a) The eigenenergies of Hq are shown in red. The black dotted lines correspond to the 
uncoupled, i.e. g = 0, case. Note that the ground state is not modified by the coupling g. The 
coupling introduces a splitting Ig^fn which is nonlinear in the phonon number n. (b) The time 
averaged squared resonator displacement as a function of the detuning of the driving frequency 
8(0 = CO — (Op for a fixed integration time of x = 100 /is. The upper (lower) blue dashed line is 
the evolution at T = OmK of the initial state | | 0) (| 10)). The strongest difference between the 
amplitudes of the two aforementioned initial states is found for a detuning of 8(0 /2k — ±0.07 MHz 
(see mark C). The upper and lower red solid lines correspond to the same initial spin states at a 
temperature of T = 30 mK. The other parameters are X/2tz = 0.04 MHz, (o p /2tz = (O q /2n — 
1.5 GHz, r = 5 • 10V 1 , and g/2n = 0.3 MHz. At 30mK the thermal energy k B T is 2.4 times 
smaller than h(O q . The letter A in both figures marks the transition from the ground state to the 
first excited state with positive detuning. The line marked B corresponds to the next transition 
between states with positive parity. A transition between the states with negative parity requires 
a negative detuning. In between A and B only two-phonon transitions can occur. The main peak 
around C at 8(0 ~ 0.07 MHz is due to several transitions with the main contribution from transtion 
between eigenstates |V4.+) — > \W5,+}- This is consistent with the magnitude of the main peak 
around n ~ X 2 /Iq ~ 4. (Inset) The squared amplitude of the CNT is shown as a function of time t 
for a detuning of the driving frequency of 8(o/2tz = 0.07 MHz. Within the first 10 ns the fast Rabi 
oscillations with the frequency of the coupling strength g are damped and following oscillations 
are caused by the driving strength A. The resonator reaches a steady state after an interaction time 

i/r. 
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a controlled state preparation is necessary. In this case ferromagnetic leads could be used as has 
already been demonstrated for nanotubes deposited on a substrate.^ 

We first solve the master equation (??) for the case of T = numerically and calculate the 
squared oscillator displacement which is proportional to the current through the charge detection 
device, in this case a QD. We assume that the oscillator frequency of the fundamental bending mode 
of the nanotube (D p /2tz = 1.5 GHz is matched by the Zeeman splitting co q /2n of the electron spin 
(the zero point motion for this case amounts to Iq = 2.5pm). At this frequency the ground state 
energy of the oscillator is 2.4 times smaller than the thermal energy at T = 30mK. Experimentally, 
frequencies of suspended nanotubes between 120 MHz and 4.2 GHz have been reported.^^H The 
coupling strength is chosen to be gj2% = 0.3 MHz which is much larger than the damping of the 
CNT, r = 0.05MHz. Together with the spontaneous relaxation of the qubit y which is negligible 
the device can be operated in the strong-coupling regime, i.e. g ^> T, y. The driving strength 
X/2tz = 0.04 MHz is chosen to be large enough to compensate for the damping but weak enough 
not to dominate over the effect of the spin-phonon coupling. 

In the inset of Figure ^p) the response of the oscillator at temperature T = is plotted as a 
function of the driving time t beginning from the initial preparation of the spin state as described 
above. The driving frequency is detuned from the resonant oscillator and qubit frequencies co p = 
co q by 8co = co — co p . The value of 8(o/2n = 0.07 MHz is found to give the largest difference in 
amplitude with respect to both initial spin orientations as we discuss below. 

A considerable difference between different initial spin states in the squared amplitude (and 
therefore in the measured current) can be observed. For the initial state |0 t), see inset of Figure 
(|2})), after the inital fast Rabi oscillations with a frequency g disappear (on a timescale rs 10ns), the 
dynamics are governed by the slower Rabi oscillations at the frequency of the drving strength A. 
The evolution of the initial state |0 \) shows only little effect of the driving. This can be explained 
by the nonlinearity introduced by the (strong) spin-phonon coupling in the Jaynes-Cummings 
Hamiltonian (??). When driving with frequencies detuned less than (y/2 — w 0.12MHz 
from resonance, the probability to leave the ground state is much less than for all other states. 
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Typically, the charge detector is too slow to follow the instantaneous motion of the resonator, 
and it will thus detect an averaged signal I °c X 2 caused by the vibrating charge. In the main part 
of Figure we plot the integrated averaged squared resonator displacement X 2 as a function of 
the detuning of the driving frequency 8(0. As an integration time we choose x = 100 /is; however, 
any integration time x of more than 30 /is was found to yield a useful signal. 

The peaks can be explained with the help of the spectrum (Figure (|2^)) of the Jaynes-Cummings 
Hamiltonian (??). Whenever the driving frequency hits a resonance, energy is absorbed by the 
oscillator which results in an increased amplitude. While only positive detuning are shown here, 
note that the result for negative detunings is exactly the same. The peak at 8(0 /2n = 0.3 MHz = 
g/2n, denoted by A, for initial state | 4- 0) for example is caused by a strong resonant coupling 
between the ground state 1 1 0) and the dressed state | Note that the smaller peaks between 

A and B are due to two-phonon processes which are less pronounced because of the weak driving. 
The main peak C centered around 0.07 MHz is due to transitions between states of equal parity, i.e. 
Wn,+) Wn+\,+) withn > 1. The line B marks the transitions — > |V2, +} which corresponds 
to a detuning of the driving frequency of 8(0 ~ 0.1 2 MHz. The main contribution to the large peak 
C comes from the transition |v / 4.+) — > |V / 5.+) which corresponds to 0.07 MHz. This is consistent 
with the peak hight of X 2 m 25 pm which in turn corresponds to a phonon number n ~ X 2 / Iq ~ 4. 
The individual transitions cannot be resolved because the lines are broadened due to the finite 
damping of the CNT Note however that a finite damping is essential to allow for non-resonant 
transitions which enable population of higher state. The large energy gap between the ground state 
and the first excited state due to the nonlinearity of the spectrum is the reason why the two spin 
states 1 1) and 1 |) can be resolved. 

In Figure ^p) the solid red lines show the integrated amplitude in the case of a finite temper- 
ature of T = 30 mK. Qualitatively, the same features are observed as for T = but in this case 
the integrated amplitude for the initial state | v(4)) ls larger because of the thermal distribution of 
oscillator in containing a small admixture of |1 1) and higher states, in addition to |0 \). 

We also note that the width of the peak is only little affected by the finite temperature and is still 
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goverend by the damping of the CNT. 

To achieve the most efficient read-out scheme we strive to increase the contrast in the integrated 
amplitude between spin up and down. A way to achieve this is to change the driving frequency 
as a function of time. This compensates for the fact that the frequencies of higher transitions are 
smaller than those of the lower transitions. From Figure (|2^) we see that the transition frequencies 
between states of same parity scales with g (\/n+T— y/n) (see above). For the sake of simplicity 
in the following we change the frequency once to demonstrate the principle, but this method could 
be optimized to further increase the contrast. To achieve the maximum contrast, i.e. the difference 
between the integrated amplitudes for the two initial states with different spin orientations, we 
switch the frequency at the point right before the squared amplitude reaches its maximum. The 
parameters used are the same as in the case of continuous driving. The time at which the driving 
frequency is switched from 1500.070MHz to 1500.036MHz is 19jUs. 

The result is shown in Figure ([3^). We clearly observe an increase in the squared amplitude of 
the resonator for the case of the initial state | y t) while the evolution of the other initial state is 
barely affected. To demonstrate the effectiveness of the driving scheme, we also plot the response 
of the system when the driving is switched off at 19/is. As expected, this causes the current 
to decay exponentially with a rate T. This pulsed driving results in lower contrast as long as 
the integration time 1 is on the same order as the time scale 1 /T given by the damping and not 
considerably longer. 

To demonstrate the feasiblility of our proposal we use a simple model of a QD capacitively 
coupled to the vibrating charge distribution on the CNT to derive an estimate for the current fluctu- 
ations. In the linear (low-bias voltage) regime, the current is given by /(Kate) = V r s rfG(V r ga te)- The 
conductance G(Vg a t e ) is assumed to be a Lorentzian as a function of the gate voltage which, due to 
capacitive coupling of QD and charged CNT, in turn depends on the squared average displacement 
X 2 of the CNT. The average displacement X itself is zero. To get an estimate of the capacitance 
we use a simple toy model in which both the QD and he charged CNT are tubes of finite length 
I and equal diameter. While this is a relatively crude model, it gives a lower limit as we show in 
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Figure 3: (a) The squared amplitude as a function of the elapsed time t is shown for the cases 
of continuous (dashed blue), pulsed (dotted yellow), and variable (solid red) driving. The same 
parameters are used as in Figure (|2]). The driving frequency ft) is changed at the time t = 19 /is 
right before the amplitude is largest, (b) The three curves corredspond to the same modes of 
driving. Each curve shows the difference between the two corresponding squared amplitudes (on 
the left v-axis) and current (on the right y-axis) integrated over an interval given by the time X in the 
x-axis, which is proportional to the average current 7+ — 7i through the QD. The variable driving 
scheme results in an increase in current by a factor of about 2. 
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the Appendix. Expanding both the capacitance as well as the current 7 up to second order in X and 
in the limit of a tube radius much smaller (1 nm) than the distance between the CNT and the QPC 
(d = 50 nm) we obtain 

W)=Io + '^ = V, d G„ (l - ^1^ ) , (6) 

where Go is the maximum conductance Go = 2e 2 /h and Co the vacuum permittivity. We assume 
AV = 0. 1 mVj^Vyd = 250 /lV,^and / = lOOnm. With these values we find a constant background 
current Iq w 20 nA. The change in current for a maximum average squared amplitude of 83 pm 
(cf. Figure ^l)) is hX 2 « — 0.83nA, which corresponds to a fluctuation in the current of about 
4.2% to be detected. In Figure ([3})) we show both the integrated difference in average squared 
amplitude as well as the corresponding current for the two different initial states for the three kinds 
of driving as. 

In conclusion, we have theoretically shown how the spin state of an electron coupled to a 
nanomechanical resonator can be read out through an adjacent charge sensor. Here, we have 
studied the case where the charge sensor is a quantum dot, but similar considerations hold for 
the case of a quantum point contact, with a slightly altered value of 1%. The read-out scheme is 
all-electrical and requires no optical access to the device nor any magnetic field gradients which is 
a promising in view of its potential for scaling to many qubits. We presented a numerical study of 
the read-out using realistic experimental parameters. We use a simplified model of the QD and the 
system geometry to demonstrate that with currently available experimental techniques spin read- 
out via the mechanical motion should be possible. For small detuning of the driving frequency 
we observe a maximum contrast between initial states with spin-up and spin-down. The contrast 
decreases with increasing temperature but is still significant at dilution refrigerators temperatures. 
We have also shown that more elaborate driving schemes in which the driving frequency is a 
function of time can lead to a large increase in contrast. Additional refinements may be used 
to further increase the amplitude of the resonator and thus the sensitiveity and speed of the spin 
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readout. 
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Appendix: Capacitive coupling of CNT and QD 

Based on the physical properties of the CNT resonator, we have calculated its spin-dependent 
mean squared resonator amplitude X 2 after external driving. We propose to read out the spin 
via the influence of the oscillating resonator on the current through a carge sensing device in the 
vicinity of the CNT. Here, we provide a simplified electrostatic model of the system which we use 
to estimate the capacitive coupling. In the following we assume the charge sensing to be performed 
by a QD (similar considerations hold for the case of a quantum point contact). 

The measured curent is given by /(Vgate) = V s dG(V gate ) with a fixed bias voltage V s d and a 
conductance G(Vg a t e )- The latter depends on the gate voltage Vg ate and it is modulated by the 
oscillating charged CNT due to the capacitive coupling to the QD. We assume the QD to be operate 
in the low-bias regime so that only a single energy level is within the bias window. In this case the 
(two-terminal) conductance G as a function of the gate voltage is given by a Lorentzian 

= W = W Ay2 +4 g -e _ Vo)2 (7) 

where the conductance at resonance is one unit of the quantum of conductance, Go = 2e 2 /h, with 
a factor of 2 accounting for the spin. The width of the Coulomb peak around the center voltage Vq 
is assumed to be AV — 0.1 mVP At 30 mK a broadening of the conductance due to temperature is 
therefore negligible. The gate voltage which depends on the displacement x in the case of a single 
charged CNT is simply given by Vg at e = — e/C(d +x). 

To find a simple analytical expressions for the the capacitance for the proposed setup, we 
use a simple toy model which captures the most essential features of our proposed setup. The 
shapes of both the charge distribution on the CNT as well as the conducting island of the QD 
are approximated as tubes of length I and diameter a. In the case of d = 50nm, a = 2nm, and 
I = lOOnm, i.e., a < d < I, the capacitance is given by C(d + x) = j ^fj^ where Co is the vacuum 
permittivity. A plate capacitor of equal area of the plates would have a smaller capacitance for all 
values of d as long as the condition d 3> a holds for the two wires. Because in the proposed setup 
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the conducting island of the QD resembles more a plate than a wire our toy model underestimates 
the actual capacitance and the result should be taken as an estimate of a lower limit. 

To get a simpler expression we want to expand Eq. (??) up to second order in x. We show that 
this is a good approximation by comparing the fluctuation in gate voltage caused by the oscillating 
CNT 5Vgate to the width of the Coulomb peak AV. As shown in Figure ([3J5), the maxium x is 
less than 10pm. For this value we find 5Vg at e = — e ( gpq^j — c{d-x) ) ~ 0.02mV which is much 
smaller than AV = 0.1 mV. Expanding Eq. (??) up to second order in x yields Eq. (??), where 
Iq is the current with the CNT at rest. For a bias voltage of V s d = 250 /iV we find Iq ~ 20 nA. 
The change in current for a maximum average squared amplitude of 83 pm (cf. Figure ([3^)) is 
I 2 X^ « -0.83 nA. In Figure we plot the difference between the currents for the states with 
initial spin orientation | and 4- where the background current Iq cancels. 
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